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We discuss results of the most accurate to-date test of the low-energy electroweak sector of 
the standard model of elementary particles. Combining previous measurements with our high- 
precision calculations we extracted the weak charge of the ^^^Cs nucleus, Qw = — 73.16(29)exp(20)th 
[S. G. Porsev, K. Beloy, and A. Derevianko, Phys. Rev. Lett. 102, 181601 (2009)]. The result is 
in perfect agreement with Qw^ predicted by the standard model, Qw^ = —73.16(3), and confirms 
energy-dependence (or running) of the electroweak interaction and places constraints on a variety of 
new physics scenarios beyond the standard model. In particular, we increase the lower limit on the 
masses of extra Z-bosons predicted by models of grand unification and string theories. This paper 
provides additional details to the Letter. We discuss large-scale calculations in the framework of the 
coupled-cluster method, including full treatment of single, double, and valence triple excitations. 
To determine the accuracy of the calculations we computed energies, electric-dipole amplitudes, and 
hyperfine-structure constants. An extensive comparison with high-accuracy experimental data was 
carried out. 

PACS numbers: 11.30.Er, 31. 15. am 



I. INTRODUCTION 

Atomic parity violation (APV) places powerful con- 
straints on new physics beyond the standard model of el- 
ementary particles [ij . The APV measurements in are in- 
terpreted in terms of the weak nuclear charge Qw, quan- 
tifying the strength of the electroweak coupling between 
atomic electrons and quarks of the nucleus. At the tree 
level the weak nuclear charge is given by a simple formula 



Qw = -N + Z {I - 4: sin^ 0w), 



(1) 



where N is the number of neutrons, Z is the nuclear 
charge, and 9w is the Weinberg angle. Since sin^ 9w is 
close to 0.23, the weak nuclear charge Qw is numerically 
close to —N. 

In Ref. [2! we reported the most accurate to-date de- 
termination of this coupling strength by combining pre- 
vious measurements [J 01 with high-precision calcula- 
tions in cesium atom. We found the result Qwi^^^Gs) — 
— 73.16(29)oxp(20)th t^] to be in a perfect agreement with 
predicted by the standard model (SM), Q^ ^ 
-73.16(3) [3. In this work we provide a detailed account 
of the calculation carried out in Ref. [1]. 

Historically, APV helped to establish the validity of 
the SM [6H8|] . While a number of APV experiments have 
been carried out the most accurate measurement 

is due to Wieman and collaborators [3I]. They deter- 
mined the ratio of i?pNc//3 = 1.5935(56) mV/cm (where 
^PNC is the parity nonconserving amplitude defined be- 
low by Eq. ^ and /3 is the vector transition polarizabil- 
ity) on the forbidden 6S1/2 781/2 transition in atomic 
Cs with an accuracy of 0.35%. This measurement does 
not directly translate into an electroweak observable of 



the same accuracy, as the interpretation of the experi- 
ment requires input from atomic theory, which links Qw 
to the signal. Qw is treated as a parameter and by com- 
bining -EpNC calculations with measurements, the value 
of Qw is extracted and can be compared with the SM 
value either revealing or constraining new physics. 

The parity nonconserving (PNC) amplitude for the 
6S1/2 7Si/2 transition in Cs may be evaluated as 



PNC 



{7Su2\Do\nPu2){nPi/2\Hw\6Su2) 



■6S, 



3i/2 

{7Si/2\Hw\nPi/2){nPi/2\E>Q\GSi/2) 



n 



Here D and Hw are electric-dipole and weak interac- 
tion operators, and Ei appearing in the denominators 
are atomic energy levels. The effective weak interac- 
tion mediated by ^-bosons averaged over quarks reads 
Uw — 75 p(r), where Gi;- is the Fermi constant, 

75 is the Dirac matrix, and p(y) is the neutron-density 
distribution. 

Interpretation of the PNC measurements requires eval- 
uating Eq. ([2|). Although the underlying theory of 
quantum electrodynamics (QED) is well established, the 
atomic many-body problem is intractable. Reaching the- 
oretical accuracy equal to or better than the experimen- 
tal accuracy of 0.35% has been a challenging task (see 
Fig. [T]). An important 1% accuracy milestone has been 
reached by the Novosibirsk [l^ and Notre Dame [l3| 
groups in the late 1980s. More recently, several groups 
have contributed to understanding sub-1% corrections, 
primarily due to the Breit (magnetic) interaction and ra- 
diative QED processes p^| - |22j (reviewed in [1^). The re- 
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suits of these calculations are summarized by the "World 
average '05" point of Fig. [1] which has a 0.5% error bar 
reflecting this progress. As of 2005, the sensitivity to new 
physics has been limited by the accuracy of solving the 
basic correlation problem. Here we report an important 
progress in solving it. 
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FIG. 1: (Color online) Progress in evaluating the PNC ampli- 
tude. Points marked Paris '86, Novosibirsk '89, Notre Dame 
'90 correspond to Refs. HI], and [13]. Point "World 

average '05" is due to efforts of several groups [l^ - [2^ on 
sub-1% Breit, QED, and neutron-skin corrections reviewed 
in Ref. [S^]. The strip corresponds to a combination of the 
Standard Model Qw with measurements [1, The edges of 
the strip correspond to ±cr of the measurement. Here we ex- 
press i?pNC in conventional units of j|e|aB {—Qw/N) x lO"'^^, 
where e is the elementary charge and qb is the Bohr radius. 
These units factor out a ratio of Qw to its approximate value, 
-N. 



We wish to evaluate accurately the sum 
Cs atom has one loosely-bound valence electron 
V — 6si/2, 6pi/2, ■•■ outside a closed-shell core 



outside a closed-shell 
Is'^ 2s'^ 2p^^2 ■■■ 5^5/2- We compute atomic wave func- 
tions, energies, and matrix elements and sum over the 
intermediate states. We solve the eigenvalue problem 
Hl'i/y) — £'^1^'^) and find atomic wave functions and en- 
ergies. Our specific scheme [1^ l25l - [27j of solving the 
atomic many- b ody problem is rooted in the coupled- 
cluster method [28|. Details will be provided in the next 
section. 

The paper is organized as follows. In Sec. [Tl] we de- 
scribe the coupled-cluster (CC) approximation including 
single, double, and valence triple excitations and present 
the results for the low- lying energy levels. In Sec. IIIII we 
evaluate the hyperfine structure constants, the matrix el- 
ements of the electric dipole moment, and the PNC am- 
plitude. An analysis of uncertainty of the PNC amplitude 
is also presented. In Sec.|lV]we extract the weak nuclear 
charge from the theoretical and experimental quantities 
and discuss implications for particle physics. If not stated 
otherwise the atomic units (h — \e\ = me = 1) are used 
throughout. 



II. COUPLED-CLUSTER APPROXIMATION 

We employ an approximation rooted in the coupled- 
cluster method [l^, |33| • The key difference compared to 
the previous CC -typ e calculations for univalent atoms 
(see, e.g., Refs. |3ll - [36| ) is our additional inclusion of 
valence triple excitations in the expansion of the clus- 
ter amplitude. We refer to this approximation as the 
coupled-cluster single, double, and valence triple (CCS- 
DvT) method. Details of our approximation may be 
found in Refs. [23l [25l - l27j . Below we briefiy recapitulate 
its main features and present numerical results. 

We choose the lowest-order Hamiltonian to include 
the relativistic kinetic energy operator of electrons and 
their interaction with the nucleus and the V'^~^ (frozen- 
core) Dirac-Hartree-Fock (DHF) potential. The single- 
particle orbitals and energies Si are found from the set 
of the frozen-core DHF equations. With the DHF single- 
particle orbitals, the second-quantized Hamiltonian reads 



H 



Ho + G 



(3) 



ijkl 



Here Hq is the one-electron lowest-order Hamiltonian, G 
is the residual Coulomb interaction, a ■ and at are the cre- 
ation and annihilation operators, and N[- • •] is the nor- 
mal product of operators with respect to the core quasi- 
vacuum state |0c}. Indices ?, j, k, and / range over all pos- 
sible single-particle orbitals, and gijki are the Coulomb 
matrix elements. 

The exact many-body state {"^v) can be represented as 
follows 



1 



1 + K 



2!' 



1*1°^), (4) 



where is the lowest-order DHF state. The cluster 

operator K is expressed in terms of connected diagrams 
of the wave operator. In our approach the operator K is 
approximated by 



K^S, + D, + S, + D,+n = 



(5) 



with the double-headed arrow representing the valence 
state. Here Sy and Dy {Sc and Dc) are the valence 
(core) singles and doubles, and T„ are the valence triples. 
Compared to the core excitations, the valence amplitudes 
involve excitations of the valence electron. A popular 
singles-doubles (SD) approximation [sij corresponds to 
neglecting non-linear terms in the expansion (jjj) and va- 
lence triples in Eq. ([S]). 

Below we present the topological structure of the equa- 
tions for valence cluster amplitudes in the CCSDvT 
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approximation. The equations in the coupled-cluster 
singles-doubles approximation (i.e., without triples) are 
presented in explicit form in Ref. [33|- A detailed tab- 
ulation of the formulas for the valence triple amplitudes 
is given in Ref. [26J. Since we do not take core triples 
into consideration, the CCSDvT equations for the core 
amplitudes Sc and Dc are the same as in {3?|. 

The total energy of the valence electron is given by the 
sum of the DHF value and the correlation energy, SEy , 

-^CCSDvT = -Sdhf + SEy . (6) 

Following the notation of Ref. [13] , we can represent the 
correlation valence energy SEy as 

6E., = 6EsD + 6Ecc + SE^t , (7) 

where correction SEsd is obtained within the SD ap- 
proach, correction SEcc comes from nonlinear CC con- 
tributions, and (5i?vT is due to valence triples. 

The topological structure for the valence singles and 
valence doubles equations may be represented as [131 

~[Ho,Sy]+6E„Sy « SD + 

Sv[Sc ® Sv] + Sv[Sc ® Sc] + 

Sy [Sc ® Dy] + Sy [Sy ® Dc] + Sy [Ty] . (8) 

-[H0,Dy]+6EyDy « SD + 

Dy [Sc «) Sy] + Dy [S c «> S c] + 

Dy [Sc ® Dy] + Dy [Sy ® D c] + Dy [S c ® D c] + 
Dy[Dc ® Dy] + Dy[Ty] . (9) 

Here [Hq, S{D)y] are commutators. Sy[Sc ($• Sy] stands 
for a contribution resulting from a product of clusters Sc 
and Sy. All other terms are defined in a similar fash- 
ion. SD terms encapsulate contributions from the SD 
approximation [3ll |. 

For valence triple amplitudes we obtain symbolically, 

- [Ho , Ty] + 6EyTy W Ty [D c] + Ty [D y] . (lO) 

Contributions Ty[Dc] and Ty[Dy] denote the effect of core 
and valence doubles on valence triples, respectively. In 
the present analysis we include only these effects, while 
omitting the effect of valence triples on valence triples 
and nonlinear CC contributions [27l |. These are higher- 
order effects that are prohibitively time-consuming for 
the 55-electron Cs atom. 

A numerical solution of the CCSDvT equations pro- 
vides us with the cluster amplitudes and correlation en- 
ergies. With the obtained wave functions for two valence 
states w and v we may evaluate various matrix elements 
(MEs), 

^wv ' — " ^ • i J 



The corresponding CCSDvT expressions are given in 
Ref. [1^. There are two important modifications com- 
pared to the earlier computations [l^: (i) explicit inclu- 
sion of valence triples in the expressions for matrix ele- 
ments and (ii) dressing of lines and vertices in expressions 
for matrix elements. The dressing mechanism [25j may be 
explained as follows: when the CC exponent is expanded 
in Eq. (jlll) . we encounter an infinite number of terms. 
The resulting series may be partially summed by consid- 
ering the topological structure of the product of cluster 
amplitudes, which may be classified using the language 
of n-body insertions. We include two types of insertions: 
particle- and hole-line insertions (line "dressing") and 
two-particle and two-hole random-phase-approximation- 
like insertions. 

Our CCSDvT code is an extension of the relativistic 
SD code [13 which employs a B-spline basis set [s^ ■ Our 
present version uses a more robust dual-kinetic-balance 
B-spline basis set [s^] as described in Ref. [i^]. This 
basis numerically approximates a complete set of single- 
particle atomic orbitals. Here, for each partial wave I we 
use 35 out of A^bas = 40 positive-energy basis functions 
generated in a cavity of radius i?cav = 75 bohr. Basis 
functions with -^max < 5 are used for single and double 
excitations. For triple excitations we employ a more lim- 
ited set of basis functions with ^max(?t,) < 4. Excitations 
from core sub-shells [4s,...,5p] are included in the calcula- 
tions of triples while excitations from sub-shells [ls,...,3c?] 
are discarded. A basis set extrapolation correction to in- 
finitely large £max, A^bas, and i?cav is added separately. 

Computations were done on a non-uniform grid of 500 
points with 15 points inside the nucleus. The nuclear 
charge distribution was approximated by p[r) — pojiX ~\- 
exp[(r — c)/a]) both when solving the DHF equations and 
evaluating weak interaction matrix elements. For ^'^'^Cs, 
c = 5.6748 fm and a = 0.52338 fm. 

Numerical results for the energies are presented in Ta- 
ble m The dominant contribution to the energies comes 
from the DHF values. Correlation corrections (^£'sd, 
SEcc, and SE^t) are dominated by the SD contribu- 
tion. We also incorporate small complementary correc- 
tions due to the Breit interaction, basis extrapolation 
((5i?extrap), and quantum-electrodynamic (QED) radia- 
tive corrections. The agreement between our ab initio 
and experimental values is at the level of 0.3% for the 65* 
state and 0.1-0.2% for all other states. 

Since the CCSDvT method is an approximation, we 
miss certain correlation effects (due to omitted quadru- 
ple and higher-rank excitations). This is the cause of 
the difference between computed and experimental en- 
ergies in Table HI To partially account for the missing 
contributions in calculations of matrix elements, we ad- 
ditionally correct the CCSDvT wave functions using a 
semi-empirical procedure suggested in Ref. [13] (see jus- 
tification in Ref. [13]) ■ In this approach, the valence sin- 
are rescaled by the ratio of experimental and 
theoretical correlation energies. A consistent definition 
of the experimental correlation energies {6Ey'^^) requires 
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TABLE I: Contributions to removal energies of 6S1/2, GP1/2, 
75i/2, and 7Pi/2 states for Cs in cm~^ in different approxi- 
mations. A comparison with experimental values is presented 
in the bottom panel. 





65'i/2 


6P1/2 


751/2 


7Pl/2 




27954 


18790 


12112 


9223 


SEsD 
SEcc 


3868 
-379 


1610 
-178 


827 
-60 


460 
-43 


SE^T 


-151 


-44 


-30 


-12 


TT^tot 

^CCSDvT 


31292 


20178 


12849 


9628 


SEBrcit"' 
SEqeo'' 
^-£'cxtrap 


2.6 
-17.6 
32.2 


-7.1 
-4.1 
15.5 


0.3 
-0.4 
7.1 


-2.5 
-0.1 
4.6 


VP ^ 
-C^cxpcrim 


31309 
31406 


20182 
20228 


12856 
12871 


9630 
9641 




"Ref. g] 


J; 'Ref. gi; '^Ref. 


Hi- 




removing the Breit and QED corrections from the exper- 
imental energy, i.e., 




E — 


^DHF — SEbicH " 


- SEqeb- 


(12) 



We will refer to results obtained using the described pro- 
cedure as "scaling" . 



III. EVALUATION OF THE PARITY 
NONCONSERVING AMPLITUDE AND 
SUPPORTING QUANTITIES 

Below we present details for the evaluation of the par- 
ity nonconserving amplitude for the 66*1/2 ~^ 75*1/2 tran- 
sition given by Eq. The CCSDvT method is an ap- 
proximation, and an important part of the entire problem 
lies with evaluating the theoretical accuracy of the com- 
puted PNC amplitude. The PNC amplitude cannot be 
directly compared to an experimental measurement. As 
seen from Eq. one needs to know the matrix elements 
of the electric-dipole operator, the energies, and the ma- 
trix elements of the weak interaction Hw- The quality 
of calculations of the dipole transition amplitudes and 
energy levels can be established by comparing them with 
experimental data, while for the matrix elements of the 
weak interaction such a direct comparison is not possible. 
Instead, we may consider the operator of the hyperfine 
interaction. Matrix elements of both the hyperfine and 
weak interaction are accumulated near the origin. There- 
fore, calculating hyperfine structure (HFS) constants for 
the low-lying states and comparing them with the ex- 
perimental data allows us to assess the quality of the 
constructed wave functions near the nucleus. 

The results of calculations of the HFS constants and 
dipole matrix elements between the low-lying states are 
presented in Tables [TT] and IIIIl For the HFS calculations 



we assumed a uniform distribution of the nuclear mag- 
netization (magnetization radius of 5.6748 fm) and used 
the nuclear ^-factor of 0.73772. We explicitly list DHF 
and SD values. The entry A(CC) indicates the change 
in the value caused by including the nonlinear terms in 
the equations for core and valence singles and doubles. 
Likewise, A(vT) and A(scaling) arise due to a subse- 
quent addition of valence triples and scaling. We also 
incorporated smaller corrections: line and vertex dress- 
ing (discussed in detail in IHj), the Breit interaction, the 
QED corrections, and the corrections due to the basis 
set extrapolation. In the lower panels of the tables we 
compare our theoretical results with the most accurate 
experimental results. We find that the discrepancies be- 
tween theoretical and experimental values for the HFS 
constants are 0.15-0.35%. For dipole matrix elements 
the theoretical values are within the error bars of the 
experiments. The uncertainty estimate of E'pnc can be 
carried out using geometric means -y/^TiSi/a^n' p, '' '' 



i/2'*" -fi/2 IJiaji 

as the relative uncertainty of this combination mimics 
the relative uncertainty of the ME of the weak inter- 
action {nSi^2\HwW Pi/2)- Deviations of these combi- 
nations from experimental data are shown in the upper 
panel of Fig. [2] We find that the standard deviation of 
our theoretical values for these combinations from the 
experimental values is 0.2%. 



75|,.2,7F|;2 



Deviation from experiment (%) 

FIG. 2: (Color online) Deviations of computed values (red 
filled circles) from experimental data (centered at zero). The 
upper panel displays combinations of magnetic hyperfine 

which mimic matrix el- 



l/2^n'-Pl/2 



structure constants yM^iS^ 

ements of the weak interaction. For these combinations, ex- 
perimental error bars are negligible compared to the theo- 
retical accuracy. The lower panel exhibits deviations of the 
computed dipole matrix elements from the most accurate ex- 
perimental results [3, [13 • 



Now we proceed to evaluating the PNC amplitude ([2]) 
by directly summing over the intermediate nPi/2 {n = 
6 — 9) states [13] ■ These states contribute 99% to the 
final £^pNC value. In Table ITVl we present the matrix ele- 
ments of the electric dipole operator and H\y, as well as 
the energy differences. Contributions from the 6 — 9Pi/2 
states to -EpNC are also listed. The matrix elements were 
computed in the CCSDvT approximation with dressing. 
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TABLE II: Magnetic-dipole hyperfine structure constants A (in MHz) for '^^^Cs. Results of calculations and comparison with 
experimental values are presented. See text for the explanation of entries. 





A(6S^ /ol 


1/2/ 


A(7S^ /■,) 


A(7P^ /o") 


DHF 


1424.9 


160.90 


391.53 


57.61 


SD 


2436.7 


310.73 


560.85 


98.34 


A(CC) 


-110.1 


-22.40 


-13.43 


-5.26 


A(vT) 


-29.6 


2.35 


-1.92 


1.18 


A(scaling) 


14.7 


2.59 


0.38 


0.41 


Complementary corrections: 










Line & vertex dressing 


-9.4 


-1.92 


-1.22 


-0.53 


Breit" 


4.9 


-0.52 


1.15 


-0.15 


QED 


-9.7' 


-0.05" 


-2.30^* 


-0.02"^ 


Basis extrapolation 


9.1 


0.71 


1.08 


0.10 


Final results 


2306.6 


291.49 


544.59 


94.07 


Experiment 


2298.16 


291.9135(15)" 


545.90(9)^ 


94.35(4)3 






291.9309(12)'' 






Difference 


0.36% 


-0.15% 


-0.24% 


-0.30% 



"Ref. [41|]; *Ref. ^S]; "Ref. ^S]; ''The QED corrections for the 75*1/2 and 7Pi/2 states were obtained by scaling those for the 

6Si/2 and 6P1/2 states; "Ref. ■'Refs. [H; "Ref. [H; ''Ref. 



TABLE HI: Reduced matrix elements of the electric dipole moment operator D (in a.u.) for ^''"'Cs. Results of calculations 
and comparisons with experimental values are presented. See text for the explanation of entries. 





|(6A/2l|Z)||65i/2)| 


|(7Pi/2||I>||6Si/2)| 


l(6Pl/2 


ll^l|7Si/2)| 


|(7Pi/2|p||7Si/2)| 


DHF 
SD 


5.2777 
4.4831 


0.3717 
0.2969 




4.4131 
4.1984 


11.009 
10.256 


A(CC) 


0.0717 


0.0058 




0.0528 


0.045 


A(vT) 


-0.0423 


-0.0302 




0.0038 


0.009 


A (scaling) 

Complementary corrections: 


-0.0123 


0.0033 




-0.0138 


-0.010 


Line & vertex dressing 


0.0036 


0.0016 




-0.0004 


0.001 


Breit" 


-0.0010 


0.0019 




0.0049 


-0.003 


QED 


0.0027''" 


-0.0028" 




-0.0043" 


0.005" 


Basis extrapolation 


0.0038 


0.0005 




0.0036 


0.005 


Final result 
Experiment 

Other results 


4.5093 
4.5097(45)'' 
4.4890(65)'' 
4.5064(47)^ 


0.2769 

0.2825(20)" 

0.2757(20)' 




4.2450 
4.233(22)^ 


10.307 
10.308(15)' 


"Ref. [iH]; ''Ref. i5i]; "Ref. 


; ''Ref. (55i; "Ref. [52] (as 


re-evaluated in Ref. ^42 
^Ref. il; 


1); ^Ref. 5] 




''Ref. 'Ref. i47i]; 
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We also used the scaling procedure. Energy differences 
were based on the experimental energies with the Breit 
and QED corrections removed. For the QSi/2, QP1/2, 
7Si/2, and 7Pi/2 states we used the following "corrected" 
energies 



Ecxp — ^-EBrcit — SEqeB- 



(13) 



The contribution of the higher-energy intermediate 8P1/2 
and 9Pi/2 states to the PNC amplitude is suppressed; for 
these states we used the full experimental energies. 

The results for the PNC amplitude are presented in 
Table |Vl Each subsequent line in the upper panel of the 
table corresponds to an increasingly more complex ap- 
proximation. We start from the SD approximation. In- 
clusion of non-linear CC terms ("CC" entry) modifies the 
SD result by almost 2%. At the next step we incorporate 
valence triples. Due to the importance of these terms we 
present a detailed breakdown of the associated effects. 
We distinguish between indirect and direct contributions 
from the valence triples. Indirect effects of triples come 
from modifying energies and single and double excita- 
tions through the Schrodinger equation. In Table |V] the 
relevant values are marked with "no vT in MEs" . The 
direct contribution arises from explicit presence of va- 
lence triples in expressions for the MEs. We also list the 
results obtained without scaling ["vT (no vT in MEs; 
pure)"] and including it ["vT (no vT in MEs; scaled)"]. 
We find that the PNC amplitude is insensitive to scaling. 
Note that a similar conclusion was drawn in Ref . [i^ . As 
the next step we replace the calculated energies in the de- 
nominators of Eq. ([2]) by the experimental energies i?corr 
as explained above. The resulting entries include the 
^'corr qualifier. We also include line and vertex dressing; 
the resulting matrix elements and detailed breakdown of 
results are listed in Table IIVI We also add contributions 
of intermediate states above 9Pi/2, including continuum, 
and contributions from core excitations. These contribu- 
tions are denoted as "n > 10" and "Core contribution" , 
respectively. Finally, the lower panel summarizes well- 
established non- Coulomb contributions such as the mag- 
netic interaction between the electrons (Breit), radiative 
(QED), and other smaller corrections. 

The accuracy of the PNC amplitude was estimated by 
comparing theoretical results for energies, dipole matrix 
elements, and magnetic hyperfine constants with high- 
precision experimental data (see Tables Ul IIIip . We find 
that the experimental energies are reproduced with an 
accuracy of 0.1-0.3%. Relevant dipole matrix elements 
are within the error bars of the experiments. Finally, 
since the hyperfine constants A are accumulated in the 
nuclear region, matrix elements of the weak interaction 
{nSi/2\Hw\n' P1/2) may be tested by forming the geo- 
metric mean ^/^^nS^J^^^nJp^ , Ref- HI]- We find that 
the standard deviation of theoretical values from experi- 
ment is 0.2%. As a test of stability of the final result, we 
also computed the main term using ab initio (i.e., with- 
out scaling) CCSDvT matrix elements and energies. The 
result, 0.8839, deviates by 0.18% from our scaled value 



of 0.8823 in Table |Vl Based on these tests, we assign 
an error of 0.2% to the main term. Finally, the "tail" 
lumps contributions of remaining excited nPx/2 states 
(including continuum) and core-excited states. The tail 
was computed using a blend of many-body approxima- 
tions and we assign a 10% uncertainty to this contribu- 
tion based on the spread of its value in different approx- 
imations. The final result includes smaller non-Coulomb 
corrections and its uncertainty was estimated by adding 
individual uncertainties in quadrature. 

Our uncertainty in £^pnc represents a two-fold im- 
provement over calculations |45| and a four-fold improve- 
ment over Ref. jT^J. Both calculations report a value of 
0.908 for the total Coulomb-correlated value, larger by 
0.9% than our 0.27%-accurate result. The reason for the 
shift in our more complete calculations is three-fold: (i) 
direct contributions of the triple excitations to matrix el- 
ements (0.3%), (ii) line-dressing of diagrams for matrix 
elements (0.3%), and (iii) consistent removal of Breit and 
QED corrections from experimental energies used in the 
scaling procedure (0.3%). Representative diagrams are 
shown in Fig. |31 



■41 



+ 



+ 




FIG. 3: Many-body diagrams responsible for the shift of 
the PNC amplitude compared to previous calculations. Top 
row; sample direct contributions of valence triples to ma- 
trix elements (wavy capped line) [2^. Bottom row: iterative 
equation for line dressing of the hole line in expressions for 
matrix elements |25| (similar equation holds for particle lines; 
exchange diagrams are not shown). 



As discussed, we make a distinction between the in- 
direct and direct contributions of to matrix elements 
[2^. Indirect effects of triples come from modifying en- 
ergies and sin gle and double excitations. In the previ- 
ous work [13, 45 1 it was approximately accounted for 
by a semiempirical scaling of single valence excitation 
(or Brueckner orbitals) to the ratio of the theoretical to 
experimental correlation energies. Direct contribu- 
tions to matrix elements, however, cannot be reproduced 
by the scaling and, moreover, require storing triples; 
due to large-memory requirements this was not done in 
Ref. [13, EB]. The size of the effect, -0.0029, is given 
by the difference between entries "vT (no vT in MEs; 



scaled, E^con)" and "vT (scaled, Scorr)" in TableEl The 
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TABLE IV: 
{nLJ; mj = 


Contribution to 
l/2\D^\n'L'J';mj, = 


-EpNC 

= 1/2). 


from intermediate states 


6-9^1/2- Dipole matrix 


elements 


are of the form 


n 


{75i/2 \D\ 
a.u. 




651/2 

10""i(-QvK/JV) a.u. 


perturbed 

a.u. 


10" 


Contribution 
""i(-gvK/iV) a.u. 


6 


1.7327 




0.05575 


-0.050949 




-1.8962 


7 


4.2071 




-0.03169 


-0.099227 




1.3435 


8 


0.3769 




-0.02118 


-0.117208 




0.0681 


9 


0.1423 




-0.01605 


-0.125993 




0.0181 


n 


{nPi/2 \D\ 6Si/2> 
a.u. 




75i/2 

10-"i(-gw/iV) a.u. 


perturbed 

a.u. 


10" 


Contribution 
■"i(-Qw/A^) a.u. 


6 


-1.8402 




-0.02697 


0.033573 




1.4783 


7 


0.1134 




0.01525 


-0.014705 




-0.1176 


8 


0.0305 




0.01024 


-0.032686 




-0.0096 


9 


0.0128 




0.00776 


-0.041471 




-0.0024 



Total 



0.8823 



line-dressing [25| was also not attempted previously. The 
line-dressing comes from resumming non-linear contribu- 
tions to wave functions, Eq. (|4]), in expressions for ma- 
trix elements. A structure of the all-order equations for 
the dressed hole lines is presented in Fig. [S] The value 
of the line-dressing correction, —0.0031, is listed in Ta- 
ble |Vl The direct contributions are most pronounced 
for the 651/2 ~ dipole amplitude, where they shift 

the value by 3%; their omission leads to a Aa deviation 
from experiment (47j . Similarly, discarding line-dressing 
shifts the theoretical values of Aqp^^^ by 0.8%. 

TABLE V: Contributions to -Epnc in different approxima- 
tions. -EpNC and A are in units of W~^^i{—Qw /N), where 
A'' = 78 is the number of neutrons in ^'^'^Cs nucleus. In the 
upper panel of the table A is the difference between the re- 
suhs given in this row and the previous row. In the lower 
panels A determines the respective contribution to iJpNC- 



Approximation 


-BpNC 


A 


"Main" term: 






SD 


0.8952 




CC 


0.8800 


-0.0152 


vT (no vT in MEs; pure) 


0.8911 


0.0111 


vT (no vT in MEs; scaled) 


0.8915 


0.0004 


vT (no vT in MEs; scaled, i?corr) 


0.8885 


-0.0030 


vT (scaled, -Ecorr) 


0.8856 


-0.0029 


Line dressing 


0.8825 


-0.0031 


Vertex dressing 


0.8823 


-0.0002 


Final main (n=6-9) 


0.8823 




"Tatl": 






n > 10 




0.0195 


Core contribution 




-0.0020 


Basis extrapolation 




-0.00006 


Total 


0.8998 




Complementary corrections: 
Breit" 




-0.0054 


QED' 




-0.0024 


Neutron skin'^ 




-0.0017 


e-e weak interaction'' 




0.0003 


Sum of corrections 




-0.0092 


Final -Epnc 


0.8906 




"Ref. gil; '■Ref. O; "Ref. 


[19!]; ''Ref. gl|,L 





IV. WEAK NUCLEAR CHARGE AND 
IMPLICATIONS FOR PARTICLE PHYSICS 

With the computed -Epnc we proceed to extracting the 
electroweak observable. The experiment [3|] determined 
the ratio of Epj^c/P = 1.5935(56) mV/cm. The most 
accurate /3 comes from a combined determination 0, |45| , 
(3 = — 26.957(51)a'^. As a result we arrive at the nuclear 
weak charge 

Qw("'Cs) = -73.16(29)exp(20)th , (14) 

where the first uncertainty is experimental and the sec- 
ond uncertainty is theoretical. Taking a weighted aver- 
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age, /? = -26.99(50)a|, of two determinations 0, liTj 
results in QwC^^Cs) = -73.25(29)cxp(20)th- Both val- 
ues are in a perfect agreement with the prediction of the 
standard model, = -73.16(3) of Ref. [5]. 

While our result is consistent with the SM, it plays 
a unique and at the same time complementary role to 
high-energy physics experiments. Our result (i) confirms 
energy-dependence (or running) of the electroweak in- 
teraction and (ii) places constraints on a variety of new 
physics scenarios beyond the standard model. 

In physics, the vacuum is never still. Each particle 
carries a cloud of continuously sprouting virtual particle- 
antiparticle pairs. The strength of the mutual interac- 
tion between two particles becomes dependent on their 
relative collision energy: at higher energies, the collision 
partners tend to penetrate deeper inside the shielding 
clouds. According to the SM, the interaction strength at 
low energies differs by about 3% from its value at 100 
GeV, see Fig. 01 For low energies, where the shielding 





1 SLAC 




_ This work 










CERN/SLAC\^ 

— 1 ^ — ' ' 1 ' — ^ 



0.01 0.1 1 10 100 



Energy scale (GeV) 

FIG. 4: (Color online) Running of the electroweak force. The 
strength of the electroweak coupling varies depending on the 
energy scale probed by an experiment. The plot shows the 
amount of variation relative to the strength at zero energies. 
The solid line is the prediction [g^] of the SM. High-energy 
experiments at CERN and SLAC have measured the strength 
of electroweak force at 91 GeV with an accuracy of ~ 0.1%. 
In 2005, a SLAC electron-scattering experiment [63| had de- 
termined the strength at 0.2 GeV with an accuracy of about 
0.5%. Our analysis of atomic parity violation probes the least 
energetic (30 MeV) electroweak interactions measured so far 
and the result is in perfect agreement with the SM. Overall, 
the predicted running of the electroweak force is confirmed 
over an energy range spanning 4 orders of magnitude. 

clouds are penetrated the least, previous analyses j22l.[23l| 
were consistent with no running. Here we improve the 
accuracy of probing these least energetic electroweak in- 
teractions. 

Compared to conventional particle-physics experi- 
ments, our result provides a reference point for the 
least energetic electroweak interactions. With our 
weak charge, we find the effective interaction strength, 
sin2 6'^f(£; ^ 0) = 0.2382(11). The resuh is in agree- 
ment with the SM value '64] of 0.2381(6). While an ear- 
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lier evidence for running of sin^ 9w has been obtained 
in the par ity violating electron scattering experiment at 
SLAC [63], the prediction of the SM was outside their 
experimental error bars. Our work provides a higher- 
confidence confirmation of the predicted running of the 
electroweak coupling at low energies. 

Notice that the relevant momentum transfer for ^•^■^Cs 
atom is just ^ 30 MeV, but the exquisite accuracy of 
the interpretation probes minute contributions of the 
sea of virtual (including so-far undiscovered) particles 
at a much higher mass scale. The new physics brought 
by the virtual sea is phenomenologically described by 
weak isospin-conserving S and isospin-breaking T pa- 
rameters Qw - Qw^ = -0.800 5' - 0.007T. At the 
Icr-level, our result implies \S\ < 0.45. The parameter 
S is important, for example, in indirectly constraining 



the mass of the Higgs particle [6^1 • Similarly, the ex- 
tra Z boson, Z'^, discussed in Ref. [1], would lead to 
a deviation ^ Qw - Qiv ~ 84(Miy/Mz^)2. Wc find 
Mz^ > 1.4TeV/c^, improving the present lower bound 
on the Z' mass from direct collider searches I66l. 
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